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Abstract. We report our first microlensing candidate NMS-El towards M31 from the data accumulated during 
the four years of Nainital Microlensing Survey. Cousin R and / band observations of ~ 13 x 13 field in the 
direction of M31 have been carried out since 1998 and data is analysed using the pixel technique jjro^osed by 
the AGAPE collaboration. NMS-El lies in the disk of M31 at a = 0''43'"33".3 and S = -h41°06'44', about 
15.5 arcmin to the South-East direction of the center of M31. The degenerate Paczyhski fit gives a half intensity 
duration of ~ 59 days. The photometric analysis of the candidate shows that it reached R ~ 20.1 mag at the 
time of maximum brightness and the colour of the source star was estimated to be {R — I)q ~ 1.1 mag. The 
microlensing candidate is blended by red variable stars; consequently the light curves do not strictly follow the 
characteristic Paczyhski shape and achromatic nature. However its long period monitoring and similar behaviour 
in R and / bands supports its microlensing nature. 
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1. Introduction 

Recent years have seen rapid advances in our understand- 
ing of the nature of dark matter in the Galaxy. Much of 
this progress has come from the studies of gravitational 
microlensing, a phenomenon of temporary amplification of 
the flux of a background star when a foreground MACHO 
passes close to our line of sight to this background star. In 
the last decade, several collaborations have focused their 
attention on Galactic bulge and Magellanic Clouds fol- 
lowing the suggestion of Paczyhski (1986) and reported 
more than one thousand microlensing events after mon- 
itoring millions of star in these directions (e.g. Alcock 
et al. 1993 for MACHO, Udalski et al. 1993 for OGLE 
and Aubourg et al. 1993 for EROS collaborations). These 
groups have ruled out most of the dark matter in the form 
of sub-stellar objects and estimated that only up to 20% 
of the halo could be made up of MACHOs in the mass 
range of O.ISMq to O.9M0 (Alcock et al. 2000., Lasserre 
et al. 2000). 

Although the conventional microlens monitoring is 
successful towards the Galactic bulge and Magellanic 
Clouds, it is fundamentally limited to only nearby regions 
where there are large number of resolved stars. In regions 
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like M31, which is at ~ 780 kpc (Stanek & Garnovich 
1998, Joshi et al. 2003), microlens monitoring is a cum- 
bersome process due to the large background of unresolved 
stars of M31 where the typical stellar density is as high 
as few hundred stars per arcsec^. To deal with this prob- 
lem, Baillon et al. (1993) proposed to monitor pixels of the 
CCD rather than the stars themselves. Any variation in 
the flux of stars due to microlensing would be determined 
in the corresponding variation in the pixel of the CCD. 
Therefore, by monitoring a pixel for a long time, one could 
detect a microlensing event. The pixel technique, which re- 
lies on the shape analysis of the light curve to detect mi- 
crolensing events, is implemented by the AGAPE (Ansari 
et al. 1997) and POINT-AGAPE collaborations (Paulin- 
Henriksson et al. 2002) and have reported more than 10 
possible microlensing events towards M31 (Calchi Novati 
et al. 2002, 2003, Pauhn-Henriksson et al. 2002, 2003). 
A similar approach but based on the image subtraction 
technique was proposed by Crotts (1992) and Tomaney & 
Crotts (1996) that has been further improved by Alard & 
Lupton (1998) and Alard (2000) by optimal image sub- 
traction using a space-varying kernel. The technique is 
widely used to detect different kinds of variable sources (cf. 
Bonanos et al. 2003; Taylor 2004). Recently the MEGA 
collaboration reported 14 microlensing candidates in their 
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Fig. 1. A 13' X 13' field of Nainital Microlensing Survey 
juxtaposed over M31. 

survey towards M31 using the image subtraction technique 
(de Jong et al. 2003). Two of them (ML-7 and ML-11) had 
already been reported by POINT- AGAPE in their survey 
using the pixel technique (PA-99-N2 and PA-00-S4 respec- 
tively). Both the groups give similar results although de- 
rived by two different techniques. 

Although pixel lensing technique allows us to detect 
microlensing events where the standard technique of mi- 
crolens detection is not applicable, it is limited to only 
bright stars that are either giants and super-giants or stars 
that are amplified enough to become detectable above the 
noise level of background flux. Since the search for gravita- 
tional microlensing events requires monitoring of millions 
of stars over a few years to yield a significant number of 
detections as well as to test their uniqueness, we started 
multi-band observations of M31 in the autumn of 1998 
under a program Nainital Microlensing Survey at the 
ARIES, Nainital, India in coUaboration with the AGAPE 
(Joshi 2004) . The details of our observations are given in 
the next section while the pixel technique and its imple- 
mentation is described in Sect. 3. The selection criteria for 
the microlensing search and our results are given in Sect. 
4 followed by the discussion and conclusions. 

2. Observations 

We started Cousins R and / band observations in the di- 
rection of M31 in 1998, with a smaU IKxlK CCD cov- 
ering an area of ~ 6' x 6' at the f/13 Cassegrain focus 
of 104-cm Sampurnanand Telescope. Since 1999, observa- 
tions have been carried out with a larger size 2Kx 2K CCD 
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Fig. 2. A histogram of the data taken in the 4 year survey. 
The total number of nights observed per observing season 
is given in the box. 

chip which covers ^ 13' x 13' encompassing the previous 
field. Each pixel covers an area of 0.37x0.37 arcsec^ of the 
sky. The target field is centered at a2ooo = 0''43™38^ and 
^2000 = +41°09'.l, at a distance of about 15 arcmin in the 
South-East direction from the center of M31. Fig. 1 shows 
the location of the target field superposed over the image 
of M31. In the 4 year observing run, we accumulated a 
total of 468 data points in R band and 383 data points in 
/ band during 141 observational nights spanning ^ 1200 
days. A histogram of the data taken in the survey is shown 
in Fig. 2. The average seeing during the observations was 
2.2 arcsec. A detailed overview of the observations is 
given in Joshi et al. (2003). 

3. Pixel method 

The basic principle of the pixel method is that if a star of 
fiux i^* is lensed in a crowded field, then by subtracting the 
original flux from the amplified flux of the star, one gets an 
increase in fiux equal to (A — 1)_F!, above the photon noise 
where A is the magnification of star. Now if the target 
star shows a variation in the flux with time, it reflects in 
the AF of the pixel. Thus by following AF with time, we 
can monitor the variation in the flux of the target star. 
Although to detect any variation in the fiux, this change 
should be significantly above the fluctuation level. 

In the point-like lens and uniform motion approxima- 
tion, the magnification factor A is related to the impact 
parameter u as given by the following Eq. 



A{t) = 



< + 



(1) 



Here io and Ie are the time 



where u^(i) 

of maximum amplification and Einstein time scale respec- 
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tively and Uq indicates the distance of closest approach in 
the units of Einstein radius. When u{t) « 1, the magni- 
fication is well approximated by 



Uo 



Then the change in flux of the target star is given by 



F A 




(2) 



t-to 



From the light curve of a resolved star, one can get three 
parameters, namely peak time (to), Einstein time (tE) and 
peak magnification (Amax) by fitting a theoretical light 
curve. However, for the unresolved stars as in the case 
of pixel lensing, there is a degeneracy between the ts, 
A-max and F^, (Gould 1996). Hence, high quality data (e.g. 
Paulin-Henriksson et al. 2003) or detection of the source 
star (e.g. Auriere et al. 2001) are required to accurately 
determine the Ie and hence the mass of the MACHO. 

3.1. Implementation of the pixel method in Nainital 
data 

Standard CCD reduction technique were used for the pre- 
processing of data using IRAF^ which includes bias sub- 
traction, flat flelding and cosmic ray removal. Wo added 
all the frames of a particular fllter taken on a single night 
and made one frame per fllter per night to increase the 
signal to noise ratio. 

3.1.1. Alignments 

To align the images, we have to flrst select a reference 
frame taken in good photometric condition with lower sky 
background and relatively small seeing. Wc chose an imago 
taken on January 05, 2000 as the reference frame. The 
average seeing during the night was 1.5 arcsec. 

Since the observations wore carried out on different 
nights, a star does not fall onto the same pixel in all the 
images. Therefore, we align all the images through rota- 
tion and shifting with respect to the reference frame. We 
achieve an accuracy better than 0.05 arcsec in alignment 
for most of the images. The photometric conditions like 
the Earth's atmospheric absorption and sky background 
are never the same for different nights. To bring all the 
images at same photometric level, we correct the images 
with respect to the reference frame using the statistical 
approach explained by Ansari et al. (1997). To overcome 
the residual gradient between two different frames due to 
reflected light, wc normalize all the images with respect 
to the median background of the reference frame. 



^ Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the National Optical Astronomy Observations 
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Fig. 3. Seeing distribution during the four year observing 
period. 

3.1.2. Seeing correction 

Most of our observations were carried out in the early 
hours of the night when sky conditions might not have 
been perfectly stable, resulting in large seeing which is 
conspicuous in Fig. 3 where we have plotted the histogram 
of seeing variation with time in the four year observing 
run. The seeing during our observations varies in the range 
~ 1 .4 to 3 .2 with an average seeing of 2 .2. The large 
seeing fluctuations during different observing periods may 
cause a false variation in the pixel light curve. To correct 
it, we make a superpixel of 7 x 7 pixeP (^ 2.5x2.5 arcsec^ 
area of sky) around each pixel so that most of the stellar 
flux falls on that superpixel^ although it only solves part 
of the seeing problem and does not correct for the loss of 
photon counts in the changing wings of the PSF due to 
highly variable seeing between the frames. We thus cor- 
rect the images by using the AGAPE 'seeing stabilization 
method' (Calchi Novati et al. 2002, Paulin-Henriksson et 
al. 2002) where we assume on the basis of empirical obser- 
vations that the difference in flux between any image and 
its median background is linearly correlated to the same 
difference for the reference image i.e. 



•e/ - 4>ref) + Pi 



(3) 



where (f)ref and (f)i are the flux of reference and i*^ images 
respectively while ^ref and are the flux of their respec- 
tive median background images. The a, and (ii are the 
seeing correction coefficients for the i*'* image. In Fig. 4, 
we show the linear correlation between {<j)ref — <}>ref) and 
i'Pcurr~4>ref) for two images, one taken in good seeing and 
other one in bad seeing conditions. The a, and Pi, which 



^ In the subsequent sections, we shall use the term 'pixel' 
instead of 'superpixel' 



4 



Y. C. Joshi et al.: First microlensing candidate towards M31 from the Nainital Microlensing Survey 



5 
Q 
< 



■9- 
< 





' / 

0= 0.66 '. , 




f= -30.5 • 

i • ^' / 

J/* •f^,/t 








/ ./• 

/ 

/ 

^ .... 1 .... 1 .... 1 .... 1 . 



-2000 

-2000 -1000 1000 2000 -2000 -1000 1000 2000 

A0,e, (ADU) 

Fig. 4. Linear correlation between the quantities A(f)tar = 

{(ptar - (ptar) and A(t)ref = {(f^ref - (j^ref)- Thc Icft and 

right graphs show the correlation for the images of good 
(1 .4) and bad (2 .9) seeing. The seeing of the reference 
image is 1 .5. Thc dashed lines is the y = x line while the 
continuous lines indicate y = a x+/3 line. 



3.2. Pixel light curves 

To derive a pixel light curve, we mask bad pixels which 
are either bright stars for which photometric analysis has 
already been carried out or have contaminated due to spill- 
over of the saturated stars. In this way, we reject about 
20% of the total pixels from the pixel analysis. 

3.2.1. Stability of light curves 

Since all the images are geometrically and photometrically 
aligned as well as seeing corrected, we expect the pixel 
light curve of non- variable objects to be stable with time. 
In Fig. 6, we illustrate the light curves of a stable pixel 
(782,579) in both R and / bands during its three year 
monitoring. The R.M.S. fluctuations along the light curve 
arc ^ 97 and ^ 108 ADUs in R and / bands respectively 
which are less than the average error bars of ~ 166 and 
~ 292 ADUs in the respective bands. This demonstrates 
that we have achieved a high level of stability in the pixel 
light curves after applying all the corrections. 
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Fig. 6. The light curves of a non-varying pixel (782,579) 
in R and / bands during the 1999-2001 observing season. 



Fig. 5. The variation of seeing correction coefficients a 
and /3 as a function of seeing. 



are calculated with a minimization procedure, shows 
a correlation with the seeing as shown in Fig. 5. After 
determining a and /? for each image, we correct them by 

inew _ 4>i~ 4>i~ A ,7 i ,\ 

<Pi = ^ <Pref (4) 

where 0"^"" is the modified flux after seeing correction. 
This correction scales all the images to the seeing of ref- 
erence image. 



3.2.2. Comparison: pixel vs photometric light curves 

To check the robustness of the pixel technique for our data, 
we compare the pixel light curve of a star with its com- 
plementary light curve derived through the photometric 
technique. However, for the comparison, we cannot take a 
very bright star, where the seeing correction used in the 
pixel technique does not work well nor a very faint star 
where the photometric error in magnitude determination 
is too large. Therefore, we chose a Cepheid oi R = 20.5 
mag having pixel coordinates (1186,595) for the compari- 
son. From the data derived with pixel technique, we obtain 
a period of 7.457±0.002 days for the Cepheid using the 
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Fig. 7. The R band phase hght curve of Cepheid 
(1186,595) derived with pixel and photometric techniques. 

Phase is plotted twice and in such a way that the minimum 
brightness falls at zero phase. The flux and the magnitude 
derived are not scaled. 



light curve. We thereafter define a bump in the light curve 
as 3 consecutive points rise above the 3 sigma level and 2 
consecutive points fall below the 3 sigma level. Since there 
is about a nine month gap between two consecutive ob- 
serving seasons, we constrain all these points to be in the 
same observing run. For each bump, we define a likelihood 
function as 

i = -/n( n ^('Z' ><^i| < <^ >i,a.)) (5) 
where 

Pm > <t>i) = r -4^exp[- ^'^~/f^^' ]# 

We chose only those pixels that have a primary bump with 
L > 300 and no secondary bump with L > 70 in the light 
curve. To exclude all those light curves that have more 
than one significant variation, we limit the ratio /"""'''"'^ 

-l^primary 

to be not greater than 0.1. We find ^ 19000 pixels that sat- 
isfy the above selection criteria. Most of them are clusters 
of pixels associated with significant physical variation. We 
chistcr these pixels and select only those pixels that have 
the highest value of L in their immediate neighbouring 
pixels. In this way we find 912 pixels characterized by the 
mono bump. To further check the fiatncss of the baseline of 
these pixels, we determine taking only those points 

which are used to determine baseline. A histogram of the 
/N for these pixels are shown in Fig. 8. To select the 
bumps with a fiat baseline, we reject all those pixels that 



period-searching program by Press & Rybici (1989). This 
value is in agreement with the period of 7.459±0.002 days 
determined by Joshi et al. (2003). This supports the ro- 
bustness of the pixel technique implemented on our data. 
For a comparison, the phase light curves obtained by these 
two different techniques are shown in Fig. 7. Nevertheless, 
it is seen that the photometric technique is better for the 
resolved/bright stars while the pixel technique is ideally 
suited in the case of unresolved/faint stars. 

4. Search for microlensing events 

The pixel method relies on the monitoring of the fiux vari- 
ation of a pixel with time. To identify candidate pixels for 
the microlensing event, we adopt the same procedure as 
defined by Calchi Novati et al. (2002). Since R band data 
have a better sampling density and photometric quality 
than that of the I band data, we used only R band data 
to select the light curves that are compatible with the mi- 
crolensing event and / band data will be used to analyse 
the microlensing candidate pixels. In the following sub- 
sections we briefly summarize the adopted procedure. 




Fig. 8. /N for the baseline flux. A vertical dashed line 
is drawn at /N = 5, above which all the pixels were 
rejected. 



4.1. Selection criteria 

Initially we define a base level which is the minimum value 
of the sliding average flux of 9 consecutive points on the 



have x^ /N > 5. This selection reduces the pixels to only 
670 that are used for the shape analysis test. We do not 
enforce any selection criteria on ti/2 as performed in some 



6 



Y. C. Joshi et al.: First microlensing candidate towards M31 from the Nainital Microlensing Survey 




other surveys (e.g. Paiilin-Hcnriksson et al. 2003, RifFeser 
et al. 2003) because of: a) observational gaps in our data 
of up to 10 days between two consecutive nights within a 
single observing season; b) the duration of the observations 
per season with an average observing density of once in 
three days; c) the total period of our monitoring (> 800 
days) that is large enough to exclude any other kind of 
variability and d) we are observing towards the outer disk 
of M31 where one generally expects long time scale events. 
Since the concentration of variable stars having longer pe- 
riod is higher as well as blending due to being greater see- 
ing, microlensing candidates suffer from inhomogeneity in 
the light curve. We therefore explicitly check each light 
curve by eye and eliminates all those pixels having a peri- 
odic variability following Paulin-Henriksson et al. (2003). 
We further fit a 5-parameter Paczyhski fit on the R band 
data and a non-linear least squares fit estimated as 



and determine the P^r degree of freedom (x^). We re- 
ject those light curves that has > 4. This left us with 
7 pixels. To further check the achromatism and symmetry 
in the selected pixels, we now use both R and / band data 
together. We find that 5 of them show significant varia- 
tion in I band while one of them shows no synchronized 
variation in R and / bands as well as a large value of its 
half intensity duration hence needs further monitoring. 
This finally gives us one pixel which we have taken as a 
genuine microlensing candidate named as NMS-El where 
NMS is an acronym of the 'Nainital Microlensing Survey'. 
A summary of the selection criteria is given in Table 1. 



Table 1 . Selection criteria to determine microlensing can- 
didates. 



Condition 


pixels left 


Masking bad pixels 
Bump detection 
Clustering of pixels 

iX^ /N) baseline < 5 

X' <4 

7-parameter fit 


80% 

~19000 

912 

670 

7 

1 



4.2. Microlensing Candidate NMS-El 

We find microlensing candidate NMS-El at the 
pixel coordinate (880,1402) corresponding to q;2000 = 
00''43™33.3^ ,52000 = -h41°06'44" with an error of ~ l" 
in both a and 5. The candidate lies at an angular dis- 
tance of 15 28 on the far side of the disk from the M31 
center and its position in our frames is shown in Fig. 9. 
The source is unresolved at its minimum brightness in our 
observations as can be seen to the left side of the figure, 
however, it is clearly visible when it reached its maximum 
brightness which is shown to the right side of the figure. 
The R and / band pixel light curves of NMS-El are shown 
in Fig. 10. A 7-parameter (i.e. full width at half maximum 
ti/2j time of maximum amplification to, impact parameter 
Wo and baseline flux and stellar flux for both the bands) 
Paczynski fit to the pixel light curves gives a half duration 
time of 59 ± 2 days with a peak magnification time 

of to = 1908 ± 1 days. 

4.3. Shape analysis and achromaticity 

A microlensing event possesses a characteristic Paczynski 
shape (Paczyhski 1986) with a flat baseline and achro- 
matic nature. The event NMS-El does not strictly follow 
symmetric shape as is evident from Fig. 10. It is seen that 
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Fig. 10. The R and / band pixel light curves of NMS-El in the upper and lower panels respectively. The continuous 
lines represent the results of the 7-parameter Paczynski fit. 



the shape of the light curve is broadened in I band in 
comparison to the R band. 

When we have seen pixels around (880,1402), we find 
a bright source at 2 .5 to 5 .5 south-west that reach 
up to -R ~ 20.6 magnitude having a mean {R — I) ~ 1.3. 
Its photometric analysis indicates that the source could 
be a red variable. We did not find any other bright star 
in our observations at the detection limit of R ^ 21.5 
mag. When we analyse the light curve of the same pixel 
in MDM data^ which is observed with the 1.3 m telescope 
and better seeing conditions, we sec a bump in the light 
curve during the 1998-1999 observing season along with 
a few overlapping points in the rising branch of the mi- 
crolensing candidate. We estimate the magnitude of the 
bump equivalent to 

R{A(t)) ~ 22.8; I{A(t)) ~ 20.7 

This star seems very close to the candidate pixel and can- 
not be resolved in our images. Although we cannot exactly 
estimate the magnitude of the star, its change in flux in 
the MDM light curve indicates that it is an extremely red 
variable. As each pixel of ~ 2 .5x2 .5 contains more than 
one thousand stars of M31, it is very likely that the nearby 

^ The pixel coordinate corresponding to NMS-El in the 
MDM target field is (1227,451). The light curve for this pixel 
was kindly provided by the AGAPE collaboration. 



bright stars and red variables contaminate the pixel flux 
which not only deforms the shape of light curve but may 
also change the achromatic nature of the event. 

To test the achromatic nature and effect of blending on 
the candidate event NMS-El, we have drawn a flux- flux 
diagram in Fig. 11. The fit (x^/ndf = 1.6) is linear for 
the upper portion of the light curve within our noise level 
but it shows a systematic increase of flux in A{(f)j) around 
the fainter side of the microlensing regime which suggests 
that the light curve of NMS-El, in general, is chromatic. 
However, keeping in mind the nearby red variable star con- 
tamination as we discussed earlier, it is not imexpccted. In 
crowded regions like M31, it is quite difficult to separate 
out the noisy contribution from nearby variable sources. 
The effect of blending becomes more severe in the I band 
where the amplitude of variability is generally higher. The 
large and variable seeing in our observations further wors- 
ens the contamination problem. Most of the pixel lens- 
ing candidates towards M31 are severely contaminated 
by the nearby red variable stars which is also reflected 
in the various candidates reported by Calchi Novati et 
al. (2002, 2003), Paulin-Henriksson et al. (2002) and de 
Jong et al. (2003). 

The physical parameters of NMS-El derived with the 
pixel technique as well its photometric analysis are sum- 
marized in Table 2. 
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Fig. 11. Fliix-fliix variation of NMS-El where the x- and 
y-axis show the increase of flux in R and / bands respec- 
tively. The points around (0,0) indicate the base of the 
light curve. A continuous line indicates the best fit linear 
regression line. 
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Fig. 12. The position of NMS-El at its maximum bright- 
ness in the H-R diagram. The solar metallicity isochrones 
by Bertelh et al. (1994) for log A = 7.5, 8.0 and 9.0 are also 
plotted. A vertical line is drawn to show the magnification 
direction. 



Table 2. Characteristics of the event NMS-El. 



a (J2000) 

5 (J2000) 

distance from M31 center 
to 

tl/2 

R 

max 

(R-I) 
xVd.f 



00^43'" 33.3" 
+41°06'44" 
15' 28" 
1908±1 
59±2 days 
20.1 mag 
1.3±0.2 mag 
1.3 



4.4. Physical interpretation of the event 

We carry out photometry of NMS-El during its maximum 
brightness phase using DAOPHOT (Stetson 1987). The 
candidate reached its peak brightness in our observations 
on December 18, 1999 having magnitudes 

i?~20.1; 7-18.8; (i? - /) - 1.3 

There can be an error of ~ 0.2 mag in the colour determi- 
nation. Assuming a distance modulus of (m— M)o = 24.49 

and extinction of Aji = 0.63, Aj = 0.47 towards our ob- 
served direction of M31 (Joshi et al. 2003), wo estimate 

Mil = -5.0; Mi = -6.1; (i?-/)o = l.l 

The position of NMS-El in the H-R diagram is shown in 
Fig. 12. We can see that NMS-El is a very bright red star 
at its peak brightness which suggests that the source must 
lie either on the giant branch or it must be a Ml-type 
main sequence (MS) star. If it is a MS star, its colour 



corresponds to an absolute magnitude of ~ 8 mag 

(Cassisi et al. 1998) hence needs a very high magnifica- 
tion (^4 

max > 10^) to reach a brightness oi Mn > —5 mag. 
Since NMS-El is a long time scale event, the latter possi- 
bility is loss likely. If it is a red giant star then we should 
be able to detect it in the Hubble Space Telescope (HST) 
images as HST is capable of detecting giant branch stars 
in M31. Unfortunately, there arc no observations available 
for this region in the HST archive. Although NMS-El is a 
long time scale event and red in colour, it is very unlikely 
that it is a Mira or any other long period variable as this 
event has been monitored for more than 800 days. This 
is one of the few pixel lensing events reported so far that 
has been tracked for such a long duration that allowed us 
to rule out any possibility of it being a variable star. 

Microlensing events that occur due to the objects in 
the halo of M31 generally have a longer time scale than 
the stars in the bulge or disk of M31 due to the observer- 
lens-source geometry (Baltz et al. 2003). The high value 
of ti/2 of NMS-El as well as its location on the far side 
of the M31 disk indicate that the event might be due to 
halo lensing. The MEGA collaboration has also reported 
two long time scale microlensing candidates in the disk of 
M31 using the differential technique i.e. ML-10 and ML-12 
with the durations of ti/2 ~ 47 and 133 days respectively 
(de Jong et al. 2003). The candidate ML-10 lies in our 
observed field at an angular distance of ~ 6 .6 from the 
position of NMS-El but shows a stable light curve in our 
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data since we started observations when this event had 
almost merged with the background noise. 

5. Discussion and conclusions 

The main difficulty with the pixel lensing search, in con- 
trast to the conventional microlensing search, is that one 
cannot usually be able to measure the Einstein time i b of 
the event from the observed pixel light curve which carries 
most of the physical information on the MACHO popula- 
tion. This is because the flux of a normally unresolved star 
in the absence of the lensing is unknown. There is however 
one event, PA-99-N1 reported by the POINT-AGAPE 
collaboration (Auriere et al. 2001, Paulin-Henrisksson et 
al. 2003), where an unlensed star has been identified in 
HST archival images. Wc lack such identification in the 
present study. Since we could not identify the source star, 
it is not possible for us to determine ts and hence the 
mass of the lens. Nevertheless, by proper modeling of the 
stellar populations which either can be source stars or 
lens objects, we can estimate various parameters like the 
timescalc of microlensing events, their rate of occurrence 
or spatial distribution (de Jong et al. 2003). However, we 
need to detect a reasonable sample of genuine microlens- 
ing events along with prior knowledge of the detection 
efficiency before making any conclusive estimates. 

We have reported our first microlensing candidate de- 
tected with the pixel technique on the far side of the galac- 
tic disk of M31. The peak brightness of the candidate 
Rmax ^ 20.1 mag coupled with its colour (J? — 7)o ~ 1.1 
and half intensity duration of ti/2 ^ 59 days lends sup- 
port to it being a red giant source while it being main 
sequence star is highly improbable. The half intensity du- 
ration along with its location in M31 suggest that the 
event NMS-El might be due to halo lensing. Since we still 
have not estimated the detection efficiency, the likelihood 
of the event and results of Monte Carlo simulations will 
be discussed in a future publication. The variable stars of 
M31 have also been the subject of numerous studies e.g. 
star formation history, stellar population etc.; we are try- 
ing to compile a catalogue of pixel variables in our target 
field. 

Here wc presented the results of the pilot observational 
campaign carried out towards M31. We have been contin- 
uing our observations of the same field using the new 2-m 
Himalayan Chandra Telescope (HCT), Hanlc, India and 
we expect to find a larger number of events which may 
enlighten us on the dark matter problem in the galactic 
halos of M31 and our own galaxy. 

Acknowledgments We would like to thank B. Paczyhski 
for his constructive comments which improved this paper. 
It is a great pleasure to thank Yannick Giraud-Heraud and 
Jean Kaplan for their kind support in initiating the project 
and helping us to get familiar with the pixel technique. DN 
is grateful to Japan Society for the Promotion of Science 
for an Invitation Fellowship. This study is a part of the 



project 2404-3 supported by the Indo-French center for 
the Promotion of Advanced Research, New Delhi. 

References 

Alard, C, & Lupton R. H., 1998, ApJ, 503, 325 
Alard, C, 2000, A&AS, 144, 363 

Alcock, C, Akerlof, C. W., Allsman, R. A., et al., 1993, Nature, 
365, 621 

Alcock, C, Allsman, R. A., Alves, D. R., et al., 2000, ApJ, 

542, 281 

Ansari, R., Auriere, M., Baillon, P., ct al., 1997, A&A, 324, 
843 

Aubourg, E., Bareyre, P., Brehin, S., et al., 1993, Nature, 365, 
623 

Auriere, M., Baillon, P., Bouquet, A., et al., 2001, ApJ, 553, 
L137 

Bonanos, A. Z., Stanek, K. Z., Sasselov, D. D., et al., 2003, AJ, 
126, 175 

Baillon, P., Bouquet, A., Kaplan, J., & Giraud-Heraud, Y., 

1993, A&A, 277, 1 

Baltz, E. A., Gyuk, G., & Crotts, A. P. S., 2003, ApJ, 582, 30 
Bertelli, G., Bressan, A., Chiosi, C., Fagotto, F., & Nasi, E., 

1994, A&AS, 106, 275 

Calchi Novati, S., lovano, G., Marino, A. A., et al., 2002, A&A, 
381, 848 

Calchi Novati, S., Jetzer, Ph., Scarpetta, G., et al., 2003, A&A, 
405, 851 

Cassisi, S., CastcUani, V., degl'Innocenti, S., & Weiss, A. 1998, 

A&AS, 129, 267 
Crotts, A. P. S., 1992, ApJ, 395, L25. 
Tomaney, A. B., & Crotts, A. P. S., 1996, AJ, 112, 2872 
de Jong, J. T. A., Kuijken, K. H., Crotts, A. P. S., et al., 2003, 

A&A, 417, 461 
Gould, A., 1996, ApJ, 470, 201 

Joshi, Y. C, Pandey, A. K., Narasimha, D., Sagar, R., & 

Giraud-Heraud, Y., 2003, A&A, 402, 113 
Joshi, Y. C, 2004, Ph.D Thesis, Kumaon University, Nainital, 

India 

Lasserre, T., Afonso, C, Albert, J. N., et al., 2000, A&A, 355, 
L39 

Paczynski, B., 1986, ApJ, 304, 1. 

Paulin-Henriksson, S., Baillon, P., Bouquet, A., et al., 2002, 
ApJ, 576, L121 

Paulin-Henriksson, S., Baillon, P., Bouquet, A., et al., 2003, 

A&A, 405, 15 
Press, W. H., & Rybici, G. B., 1989, ApJ, 338, 277 
Riffeser, A., Fliri, J, Bender, R., et al., 2003 ApJ, 599, L17. 
Stetson, B., 1987, PASP, 99, 191. 

Stanek, K. Z., & Garnavich, O. M., 1998, ApJ, 503, L131 
Taylor, S. F., 2004, PASP, 116, 1126 

Udalski, A., Szymanski, M., Kaluzny, J., et al., 1993, Acta 
Astron., 43, 289 



